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A  RESUME  OF  STOCHASTIC,  TIME- VARYING,  LINEAR  SYSTEM  THEORY 
WITH  APPLICATION  TO  ACTIVE-SONAR  SIGNAL-PROCESSING  PROBLEMS 

by 

Lewis  Meier 


ABSTRACT 


This  report  summarizes  linear  time-varying  systems  theory  and  its 
application  to  active  sonar  signal  processing  when  certain  of  the  para¬ 
meters  are  random  variables.  To  keep  the  analysis  within  reasonable  bounds 
the  random  parameters  are  described  by  their  first-  and  second-order 
moments  only.  The  report  starts  by  outlining  the  theory  and  follows  by 
describing  two  possible  applications:  one  for  the  derivation  of  the 
scattering  function  of  a  moving,  turning  line  target;  the  other  in  the  use 
of  pulse  trains  to  measure  scattering  functions. 
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INTRODUCTION 


This  report  is  the  sequel  to  an  ec'lier  report  <1>  on  linear,  time-varying 
systems  theory  and  its  application  o  active-sonar  signal  processing.  The 
major  concern  in  that  report  was  with  deterministic  systems,  while  the  main 
concern  in  this  report  is  with  stochastic  systems.  The  main  results  of  the 
earlier  report  may  be  summarized  as  follows:  a  linear,  time-varying  system 
may  be  viewed  as  spreading  a  signal  in  time  and  either  frequency  or  time 
stretch  (doppler);  therefore  it  is  convenient  to  represent  such  systems  by 
a  function  of  time  shift  x  and  frequency  shift  v  or  time  stretch  a 

known  as  the  spreading  function  S  or  and  to  represent  signals  by 

another  function  of  x  and  v  or  a  known  as  the  cross-ambiguity 

function  <}>  or  The  cross-ambiguity  function  of  the  output  of  such  a 

system  is  a  modified  double  convolution  in  x  and  v  or  a  between  its 
spreading  function  and  the  cross-ambiguity  function  of  the  input;  further¬ 
more  the  spreading  function  of  the  concatenation  of  two  such  systems  is  the 
same  modified  double  convolution  of  their  spreading  functions. 

In  this  report  we  are  concerned  with  linear,  time-varying  systems  for  which 

S(x  ,v)  or  S^(x,a)  are  random  variables;  such  systems  are  then  des¬ 
cribed  by  their  statistical  properties.  The  general  range  of  statistical 
behaviour  possible  is  far  too  vast  for  this  brief  report;  therefore  we 
must  narrow  our  concern  by  the  use  of  assumptions  that  vastly  simplify  the 
analysis  while  only  marginally  restricting  its  application.  It  is  quite 
reasonable  to  assume  that  the  systems  in  which  we  are  interested  are  sta¬ 
tistically  independent  of  one  another  and  their  inputs.  With  this  assump¬ 
tion  it  is  also  quite  reasonable  to  assume  that  the  means  E[S(x,v)]  or 

E[S°(x,a)]  are  zero  (where  E[  ]  signifies  expected  value  of)  since  the 
effect  of  non-zero  means  may  be  analyzed  separately  with  the  deterministic 
theory  of  <1>.  Since  we  are  dealing  with  linear  systems,  it  is  also 
reasonable  to  restrict  our  concern  to  second  order  statistics  —  the 

variance  of  S(x,v)  and  SD(x,a).  Finally  we  assume  the  values  of  S(x,y) 

or  S^(x,a)  are  uncorrelated  for  differing  x  and  v  or  a  ;  hence  the 

second-order  statistics  of  S  of  can  be  represented  by  the  functions 

$(x,v)  or  $D(x,a)  ,  known  as  the  scattering  function,  which  is  the 

analogue  of  the  power  density  of  a  white-noise  process: 

E[S*(x‘  ,v' )  S(x,v)]  =  $(t,v)6(x'-x)  6(v’-v)  , 

E[SD*(x',a')  S(x,a)]=  $D(x  ,a)  6(x'-x)  6(a'-a)  ,  (Eq.  1) 


(where  *  represents  complex  conjugate  and  6(  )  is  the  Dirac  delta 

function). 
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The  transmitted  signal  is  assumed  to  be  deterministic;  hence  its  ambiguity 
function  is  real.  On  the  other  hand,  once  the  signal  has  passed  through  a 
stochastic  system  of  the  type  discussed  in  the  last  paragraph  it  becomes 
stochastic  with  zero  mean;  therefore  in  our  stochastic  theory  signals  will 

be  represented  by  power  cross-ambiguity  functions  E[|<|»|2]  or  E[|<t^|2]  . 

The  fundamental  result  of  import  is  that  the  power  cross-ambiguity  function 
of  the  output  of  a  system  describable  by  a  scattering  function  is  a 
modified  convolution  in  t  and  v  or  a  between  its  scattering  function 
and  the  input  power  cross-ambiguity  function;  furthermore  the  scattering 
function  of  the  concatenation  of  two  systems  describable  by  scattering 
functions  is  the  same  modified  convolution  in  x  and  v  or  a  between 
their  scattering  functions.  In  fact,  the  required  modified  double  convolu¬ 
tion  in  x  and  v  is  just  a  double  convolution,  while  the  required 
modified  double  convolution  in  t  and  a  is  the  same  modified  double 
convolution  required  by  the  deterministic  theory  and  discussed  in  detail  in 
<1>.  These  results  are  illustrated  in  Figs.  1  and  2;  the  frequency  shift 
version  of  the  results  is  well  known  <2>,  but  the  time  stretch  (doppler) 
version  is  believed  to  be  novel. 

The  fact  that  in  the  frequency-shift  version  of  the  theory  an  ordinary 
double  convolution  is  involved  should  not  be  a  big  surprise.  The  assump¬ 
tion  that  Eq.  1  holds  for  S(x,v)  implies  that  E[h(t',x')  h(t,x)]  is 
non-zero  only  for  x'  =  x  and  depends  only  on  t1 -t  ;  that  is,  that  the 
time-varying  impulse  response  h(t,x)  consists  of  a  statistically  wide 
sense  stationary  set  of  uncorrelated  scatterers.  Often  this  assumption  is 
abbreviated  WSSUS  for  wide  sense  stationary  uncorrelated  scatterers.  The 
implication  of  this  result  is  that,  unlike  the  deterministic  case  and  the 
time-stretch  version  of  the  stochastic  case,  the  order  in  which  systems  are 
concatenated  is  unimportant.  These  results  are  exactly  analogous  to  the 
results  of  time-invariant  linear  system  theory  in  which  ordinary  convolu¬ 
tions  are  used  and  the  order  of  concatenation  is  unimportant.  Unfortunate¬ 
ly,  while  the  frequency-shift  scattering  function  is  a  good  model  of  the 
medium,  to  correctly  represent  a  thin,  rigid,  moving,  turning  target  (such 
as  a  submarine)  the  time-stretch  scattering  function  is  required  (and  even 
this  is  an  approximation).  For  unity  BT  signals  such  as  CW  signals,  such  a 
target  may  be  represented  approximately  by  a  frequency-shift  scattering 
function  derived  in  an  obvious  way  from  the  time  stretch  scattering 
function,  but  for  high  BT  signals  such  as  linear  FM  signals  the  time- 
stretch  version  must  be  used. 

In  many  cases  it  is  desirable  to  try  to  measure  scattering  functions.  For 
example,  knowing  the  scattering  function  of  a  medium  tells  us  the  statis¬ 
tics  of  how  it  spreads  in  time  and  frequency  a  signal  transmitted  through 
it,  and  knowing  the  scattering  function  of  a  target  allows  us  to  measure 
interesting  target  properties  <3>.  Thus  the  design  of  a  suitable  signal 
for  measuring  scattering  functions  is  of  considerable  interest  and  has  been 
the  subject  of  considerable  studies,  for  example  by  Rihaczek  <4>.  The  rub 
is  that  the  volume  under  the  power  ambiguity  function  |y|2  is  constrained 
to  a  given  value,  so  that  it  is  impossible  to  determine  the  scattering 
function  exactly.  Costa  and  Hug  <5>  suggested  the  use  of  pulse  trains  to 
measure  the  under-spread  scattering  functions  (those  whose  extent  in  time 
and  frequency  has  a  product  less  than  unity).  A  pulse  train  has  a  sharp 
spike  at  the  origin,  with  the  remainder  of  its  volume  spead  outside  a 
region  whose  area  is  unity.  For  an  evenly  spaced  train  of  identical  pulses 
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this  volume  consists  of  respective  spikes  in  both  the  time  and  frequency 
shift  directions,  but  these  may  be  modified  by  any  number  of  tricks  such  as 
are  given  in  <4>.  In  their  work  Hug  and  Costa  make  implicit  use  o.  an 
approximate  method  of  realizing  the  bank  of  matched  filters  required  to 
determine  the  ordinary  output  power  cross- ambiguity  function  that  can  most 
simply  be  described  as  the  temporal  analogue  of  narrowband  (phase-shift) 
beamforming.  In  this  report  a  similar  technique  is  described  for  realizing 
the  bank  of  matched  filters  required  to  determine  the  time-stretch  output 
power  cross-ambiguity  function  that  is  the  temporal  analogue  of  broadband 
(time-shift)  beamforming. 

This  report  is  divided  into  two  parts  covering  theory  and  applications 
respectively.  The  theoretical  part  (Ch.  1)  is  largely  a  detailed  presenta¬ 
tion  of  the  basic  relationship  given  in  Figs.  1  and  2,  while  the  part  on 
application  (Ch.  2)  derives  the  time-shift  scattering  function  of  a  moving, 
turning,  line  target  and  investigates  the  use  of  coherent  pulse  trains  in 
measuring  target  scattering  functions.  Detailed  derivations  are  relegated 
to  the  appendices. 


FIG.  1  BASIC  INPUT-OUTPUT  RELATIONSHIP 
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1  THEORY 

This  chapter  consists  of  two  parts:  discussion  on  scattering  functions 
(Sect.  1.1)  and  discussion  of  the  fundamental  relationships  (Sect.  1.2). 


1. 1  Scattering  functions 

The  scattering  function  of  a  stochastic,  linear,  time-varying  system  was 
given  in  Eq.  1  in  terms  of  the  spreading  function  of  such  a  system,  where 
the  spreading  function  is  defined  in  detail  in  <1>.  Not  all  stochastic, 
linear,  time-varying  systems  have  scattering  functions,  of  course,  since  a 
number  of  assumptions,  which  were  enumerated  in  the  introduction,  have  to 
be  satisfied;  furthermore,  if  a  system  has  e  frequency- shift  scattering 
function  it  will  not  have  a  time- stretch  scattering  function  and  vice- 
versa.  Nevertheless,  many  important  systems  come  close  enough  to  obeying 
the  required  assumptions  that  they  can  be  represented  by  scattering 
functions  and  in  some  situations  a  system  with  a  time-stretch  scattering 
function  may  be  represented  by  a  frequency-shift,  scattering  function. 

Recall  from  <1>  that  the  impulse  response  of  a  linear,  time- varying  system 
is  related  to  its  spreading  function  via 


00 

h(t,x)  =  J 


e2njv(t-t) 


-00 


5(t,v )  dv 


(Eq.  2) 


therefore  if  the  system  has  a  scattering  function 


E[h*(t' , t  1  )h(t,t)]  =  J  S  e27Tj[v(t"T)'v'(t'"xl)]ElS*(T'  ,v '  )S(x ,v)]dv  dv' 


=  X  e2njv(t"t,)$(x,v)dv  6(x'-t)  (Eq.  3) 


Thus  h(t,t)  consists  of  uncorrelated  scatters  whose  statistics  are  time 
invariant.  Note  the  presence  in  Eq.  3  of  the  inverse  Fourier  transform  of 
the  scattering  function.  Just  as  in  <1>  there  was  no  simple  relationship 

between  S^(x,v)  and  h(t,x)  ,  there  is  no  simple  relationship  between  the 

statistics  of  h(t,  x)  and  $^(r,a)  for  systems  having  a  time-stretch 
scattering  function. 
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1. 2  The  fundamental  relationships 

Consider  the  two  operators  (x)  and  ®  ' 

t  ,v  i  ,a 


defined  by 


CO  00 

b  (x)  a(t,v)  =  j'J  a(x-t v-  v)‘  )  b(T*  ,v  1 )  dV  dx 1  (Eq.  4a) 

I  ,V  -<X>  -CO 

00  00 

bD  ®  aD(x  ,a)  =  J  J  aD[a'  (x  -  x') ,  ^,]  bD(x  1  ,a' )  do'  dx'  (Eq.  4b) 
x  ,a  -  oo  ° 


The  former  is  nothing  more  than  the  well-known  double  convolution  operator, 
which  is  both  communicative  and  associative.  The  latter  is  the  modified 
double-convolution  operator  studied  in  detail  in  <1>,  where  it  was  shown  to 
be  associative  but  not  communicative.  Use  of  these  definitions,  those  of 
Eq.  1  and  the  fundamental  input-output  and  concatenation  relationships  for 
spreading  and  cross-ambiguity  function  given  in  <1>  readily  yields  the 
fundamental  relationships  described  in  Figs.  1  and  2  -  details  are 
contained  in  Appendix  A. 
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2  APPLICATIONS 

Two  applications  of  the  theory  of  Ch.  1  are  presented  here:  derivation  of 
the  scattering  function  of  a  moving,,  turning,  line  target,  and  the  use  of 
pulse  trains  to  measure  scattering  functions.  Discussion  of  the  latter 
topic  is  broken  in  turn  into  two  parts:  derivation  of  its  power  ambiguity 
function,  and  presentation  of  an  approximate  method  of  simply  realizing  the 
matched  filter  bank  when  a  pulse  train  is  used. 


2. 1  The  doppler  scattering  function  of  a  moving,  turning,  line  target 


From  Eq.  22  of  <1>  the  doppler  spreading  function  of  a  moving,  turning 
target  is 


SD(x,a)  = 


c 

( 1  -  a )  6 


■f 


■CT+H 


ot-1 

a 


(1-a)  6 


-  i/c 


^?)]  • 


(Eq.  5) 


where  f(a)  is  the  reflectivity  of  the  point  along  the  target  located  at 
distance  a  from  the  centre.  The  remaining  quantities  in  Eq.  5  are 
defined  in  <1>  and  since  they  are  not  of  direct  concern  their  definition  is 
not  repeated  here.  Now  we  make  the  assumption  that 


E[f  (a')f(a))  =  Rf(  a)  6(a'~a) 


(Eq.  6) 


For  our  theory  to  be  valid  it  is  not  necessary  that  Eq.  6  hold  exactly,  but 
only  that  the  reflectivity  of  points  on  the  target  be  uncorrelated  at 
distances  greater  than  the  distance  resolution  of  the  transmitted  signal  -- 
a  condition  that  is  usually  met  in  practice. 


Now  from  Eqs.  5  and  6  we  have 


CD*  ,  ,._D, 

b  ( t  ,ot  )S  (t 


.«)13 


n%5 


„  r  ,  Ya'-l  a-1  Vi 

M  (T -a)-i-  )  T  -T+h^  •  —  )\ 

•  «[•-’  4  -  &  (- ■  k*  ^ ^)] 

irr-’-k0'1  + 1  -  ^?(T'  I +  ^  ^)] 


!f(: 


-CTt£--_— 


•  6(a'-a)  , 


(Eq.  7) 
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making  use  of  identities  implied  by  the  delta  function.  Obviously  from 
Eq.  1  and  this  result,  the  doppler  scattering  function  is: 


-CT  +  Z 


g-1 

a 


(!-*)« 


Or} 

a 


* 


(Eq.  8) 


which  has  the  same  form  as  that  of  the  spreading  function  as  given  in 
Eq.  5. 


2.2  The  power  ambiguity  function  for  a  pulse  train 

In  this  and  the  following  section  we  are  concerned  with  an  evenly  spaced 
train  of  identical  pulses: 


XPT  (t)  =  ±  ZXp[t-  (i  -^)  AT]  (Eq.  9) 

where  the  pulse  XD(t)  is  assumed  to  be  normalized  to  have  unity  total 
energy  (i.e.  the  integral  of  |xp|2  over  t  is  unity). 

Rihaczek  <4>  has  already  determined  the  ordinary  ambiguity  function 
Tpj(t,v)  for  such  a  signal  to  be 


n-1 

Ypr(x.v)  =  Z  ^(x-iAT.v)  ,  (Eq.  10) 

i=-n+l 


where 


Y-j  (x ,  v) 


e~injvAT 


sin[nv(n-  lil  AT] 
n  sin(n  v  AT) 


Yp(x,v) 


(Eq.  11) 


and  y  (t,v)  is  the  ordinary  ambiguity  function  of  Xn(t). 

P  P 

If  we  assume  —  as  is  typically  the  case  —  that  the  spread  of  yp  in  x 

is  small  compared  with  At  ,  the  ordinary  power  ambiguity  function 
| Ypy (x , v) | 2  for  the  pulse  train  is 

n-1 

|Ypt(x,v)  |2=  -1  Iy^x.v)!2  , 


(Eq.  12) 
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where 

I  Vi  Cx  ,v)|2 


sin[nv(n-  [i|  )AT]  2 
n  sin(nvAT) 


!  Vp(t,v)|2 


(Eq.  13) 


Except  for  |  i  |  =  n,  the  first  factor  on  the  right  side  of  Eq.  13  is  a 
function  of  v  that  is  a  series  of  spikes  that  are  (n- |i  |)/n  tall, 
roughly  l/[(n-  i  )AT]  wide,  and  separated  by  1/AT  .  Since  the  width 
of  ! Vp  1 2  is  small  compared  with  AT,  I  Ypjt  2  is  thus  a  pattern  of  spikes 

in  (x,v)  space.  This  pattern  is  discussed  and  illustrated  in  detail  in 
<4>. 

Now  consider  the  doppler  power  ambiguity  function  of  the  pulse  train,  in 
Appendix  B  it  is  shown,  by  substitution  of  Eq.  9  into  the  definition  of  the 
doppler  ambiguity  function  given  in  <1>,  that 


D  r\ 

Y  (x,a)  =  I  VU(t-i  AT,  a) 

PT  i=-n+l  i 


(Eq.  14) 


where 


D,  .  1  "'l1 I  D  a  -1  n+1 

Yi<T'“)=K  ^  Yp  T  *  a  T 


k  + 


l-i 


AT,  a 


(Eq.  15) 


aid  Vp(T,a)  is  the  doppler  ambiguity  function  Xp(t)  .  The  analogy  between 
Eqs.  10  and  14  is  immediate;  Eqs.  11  and  15  are  best  compared  through 
Fourier  transforms,  which  are  obvious  for  y^  and  derived  in  Appendix  B 
for  y i  : 


r1(f,v)  =  e 


rlj(ffO)  =  e 


•injvAT  ( sinfnv(n-|i|)ATl\  r  . 
V  n  sin(TtvAT)  )  V  ‘  J 


-i»j  2llfAT/sin  [*  2i1f(r'-ht)AT] 
n  sin  f  AT^ 


Tp(f,a)  , 


(Eq.  16) 


(Eq.  17) 


where  r  and  r  are  the  Fourier  transforms  with  respect  to  x  of  y 
P  P  p 


and  y  respectively. 
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Again  we  assume  that  the  spread  of  y  intis  small  compared  with  AT 

r 

In  this  case,  the  doppler  power  ambiguity  function  |vpj(t,a)|  2  for  the 
pulse  train  is 


(t,v  ) 


n-1  D 

*  ]X  (t,v)|  2 

i=-n+l  i 


(Eq.  18) 


for  |a-l  |  /a  <  1/n  .  This  constraint,  which  follows  direct  from  Eq.  15 
and  is  derived  in  Appendix  B,  insures  that  the  y?  do  not  overlap  significant¬ 
ly.  The  analogue  of  Eq.  13  is  not  given  since  the  magnitude  squared  of 
Eq.  15  does  not  take  a  simple  form. 

Instead  we  determine  the  total  energy  in  y^.  (t,a)  for  a  given  a  by  integra- 
ting  |y-(x,a)|  over^  t  or  equivalently  from  Parseval's  relationship  by 
integrating  |  f.  (f,a)|  over  f  .  If  a  similar  procedure  is  carried  out  for 
Eq.  16  the  result  is  (sin[nv(n-| i | )  AT]/[n  sin(nvAT)])2  multiplied  by  the 
energy  in  y  ''t.v)  for  a  given  v  .  For  a  bandpass  pulse,  rpD(f,a)  will  be 
centred  more  or  less  around  f  =  fQ,  the  carrier  frequency;  therefore 
the  result  of  integrating  |  r 9 ( f  ,a ) |  is 

{ s i n [ 7t  ^  fG(n-|  it  )At]/[n  sin(n  ^  f0AT)]}2 

smoothed  over  a  ,  multiplied  by  the  energy  in  y^(t,a)  for  a  given  a  . 

D  2  ^ 

Thus  |Tpj|  is  a  pattern  of  spikes  in  (t,v;  space  analogous  to  the  pattern 

of  spikes  in  (t,v)  space  of  | Vp j 1 2  ;  however,  since  the  smoothing  effect 

just  mentioned  increases  linearly  with  a  ,  the  spikes  become  lower  and 
broader  with  increasing  a  . 

To  illustrate  the  ideas  of  the  previous  paragraph,  an  example  will  be 
presented.  Let  Xp(t)  be  an  FM  pulse  with  centre  frequency  fc  and 

bandwidth  B  .  From  <1> 

I  FpM(f,a)|2  ft  |  [  f  -  fo  ■  (B  -|  c-il  fo)/2]  ,  (Eq.  19) 


where 


•■T>  ={o 


1 1  5  T 
otherwi se 
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therefore 


-co 


(B-|a-l|f0)/2 


(Eq.  20) 


where 


4 


_  g-1 
a 


fc  AT 


The  integral  of  Eq.  20  is  readily  performed  since  the  integrand  can  be 
written  as  a  sum  of  cosine  functions.  On  the  other  hand  for  the  ordinary 
power  ambiguity  function  we  have 


(Eq.  21) 


where  4=VAT  and  p  is  as  before.  Graphs  of  these  two  functions  versus 
4  are  given  in  Fig.  3  for  p  =  0.06,  |a-l|,  and  v  /fQ  «  p  ,  i  =  0  , 
and  n  =  2,  3  and  5.  These  graphs  bear  out  the  conclusions  of  the  previous 
paragraph. 

From  the  foregoing  it  is  apparent  that  evenly  spaced  trains  of  identical 
pulses  have  ambiguity  functions  with  a  central  spike  surrounded  by  an 
essentially  clear  region  of  unity  area;  hence  by  appropriate  choice  c' 
signal  parameters  it  is  possible  to  estimate  a  scattering  function  known  to 
be  contained  in  region  of  unity  area. 


13 
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FIG.  3  NORMALIZED  ENERGY  IN  AMBIGUITY  FUNCTION  FOR  i  =  0 
a)  n  =  2  b)  n  =  3  c)  n  s  5 
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The  cost  of  such  capability  is  ambiguity  about  the  absolute  position  of  the 
scattering  function.  In  many  cases  -  for  example  determination  of  medium 
scattering  functions  -  this  ambiguity  is  unimportant  since  it  may  be 
resolved  by  other  information.  It  is  possible,  as  Rihaczek  <4>  discusses 
at  length,  to  modify  an  evenly  spaced  train  of  pulses  in  such  a  manner  that 
the  non-central  spikes  of  the  ambiguity  function  are  lowered  and  spread 
out.  Without  going  into  details,  tne  possibilities  are  as  follows:  By 
shifting  the  phase  of  the  individual  pulsesDor,  better  yet,  making  the 
spacing  between  pulses  non-uniform  or  y.  may  be  lowered  and  spread 

in  t  -  thus  alleviating  the  ambiguity  in  x  .  By,  in  addition, making  the 
pulses  non-identical  -  for  example  stepping  them  in  frequency  -  it  is 
possible  to  generate  a  signal  with  a  "thumb  tack"  ambiguity  function  (i.e. 
a  central  spike  surrounded  by  a  more  or  less  uniform  plateau).  The  cost  of 
the  latter  s':ep  is  that  the  simplified  processing  for  pulse  trains  of 
identical  pulses  to  be  described  next  cannot  be  used. 


2 . 3  An  approximation  to  the  matched  filter  bank  for  a  pulse  train 

For  ease  of  presentation  a  uniformly  spaced  train  of  identical  pulses  is 
assumed;  however,  the  same  analysis  applies  in  an  obvious  manner  if  the 
pulses  are  shifted  in  time  or  phase.  Furthermore,  it  is  assumed  that  the 
individual  pulses  are  FM  pulses;  however,  the  same  analysis  again  applies 
in  an  oL  /iois  manner  to  any  wide  band  pulses. 

Let  y(t)  oe  the  signal  to  be  passed  through  the  bank  of  filters  matched 
to  Xpj(t)  and  frequency-shifted  versions  of  xpy(t).  From  Eq.  9  and  the 

approximate  identity 


foV 


2n  j  vt. 

e  x 


-2nj 


FM 


(t) 


x  +  i?) 

FM  K 


(Eq.  22) 


where  f0  is  the  centre  frequency  and  k  the  frequency  slope  of  Xp^Ct)  , 

it  is  shown  in  Appendix  C  by  substitution  that  the  output  of  the  bank  of 
matched  filters  is 


4>pt(t  ,v)  =  J  e 


f  v 

2"J-r 


-2njv(t-x)  * 


XpT(t-x)  y(t)  dt 


n  -2irjv(i-  ^>r-)AT 
E  e  c 


/I T  i  =1 


^FM^t"  if  +  (1-5F-)aT»  o] 


(En.  23) 
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where 


4»fm(t,o)  =  J  XpM(t-T)  y(t)  dt  (Eq.  24) 

—00 

Now  consider  a  bank  of  filters  matched  to  Xpy(t)  and  doppler-shi fted 
versions  of  xpT(t).  From  Eq.  9  and  the  approximate  identity 

XFM(at)  *  e  Xp^[t  +  (or  1)  fc/k]  ,  (Eq.  25) 

it  is  shown  in  Appendix  C  by  substitution  that  the  output  of  the  bank  of 
matched  filters  is 


4>pT(t,a)  =  J  X*  [a(t-x)]y(t)  dt 

-00  py 


(o-l  )f02 
2nj  - t - 


(Eq.  26) 


The  quantity  <Jifm(t,o)  is  nothing  more  than  the  output  for  a  filter 
matched  to  the  FM  pulse  Xp^(t);  therefore  the  first  step  in  the  approx¬ 
imate  realization  of  either  the  frequency-shift  or  doppler-shift  bank  of 
matched  filters  is  to  pass  the  signal  through  a  filter  matched  to  the 
transmitted  pulse.  For  a  given  x  ,  Eq.  23  for  the  frequency-shift  bank  of 
matched  filters  takes  -  except  for  a  time  and  phase  shift  -  the  form  of  a 
discrete  Fourier  transform  of  the  sequence 

<t>FMU+(i  -  ^  )  AT,  o]  ,  1  =  1,  ....  n 


*FM  'T 


(ct.  -VL  fiL  ♦  \  ( i  "  AT,  o] 


of  Outputs  from  the  initial  matched  filter.  Similarly  for  a  given  t  , 
Eq.  26  for  the  doppler-shift  bank  of  matched  filters  takes  -  except  for  a 
time  and  phase  shift  -  a  form  similar  to  the  discrete  Fourier  transform  of 
VFM  except  that  time  stretch  by  1/a  replaces  frequency  shift  by  -v  (as  might 
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1 

be  expected).  In  fact  if  "Cl  -  -jp)AT,  o]  is  identified  as  the 

matched  filter  output  for  a  hydrophone  located  at  point  (i  -  ^—)  Ax  along 

a  linear  ar'ay,  Eqs.  23  and  26,  with  At  replaced  by  Ax,  take  the  form  of 
a  narrow-band  beam  former  and  a  broadband  beamformer  respectively. 

From  this  correspondence  it  is  clear  that  Fig.  3  gives  narrow-band  and 
broadband  beam  patterns  for  2,  3  and  5  element  arrays  where  the  abscissa  is 
the  cosine  of  the  bearing  angle  relative  to  the  antenna  suitably 
normalized. 
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CONCLUSIONS 


Under  appropriate  conditions  random  linear  time- varying  systems  may  be 
represented  by  a  function  of  time  shift  and  frequency  shift  or  time  stretch 
known  as  the  scattering  function  and  the  ouptut  of  such  systems  by  a  second 
function  of  the  same  variables  known  as  the  power  cross- ambiguity  function. 
The  input-output  and  concatenation  relationships  for  such  systems  are 
double  convolutions  in  the  frequency-shift  version  and  modified  double 
convolutions  in  the  time- stretch  (doppler)  version.  A  moving,  turning  line 
target  has  a  doppler  scattering  function  that  is  non- zero  only  along  a  line 
in  time-shift,  time-stretch  space.  A  very  useful  signal  to  estimate 
scattering  functions  such  as  those  of  a  medium  or  a  target  is  the  pulse 
train,  because  its  ambiguity  function  has  a  spike  at  the  origin  surrounded 
by  a  clear  region  of  unit  area.  The  matched  filter  bank  for  such  a  signal 
when  the  individual  pulses  are  identical  and  wide  band  consists  of  a  filter 
matched  to  the  individual  pulse  followed  by  the  temporal  analogue  of  a 
narrow-band  or  broadband  beam  former. 
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APPENDIX  A 


DERIVATION  OF  THE  FUNDAMENTAL  RELATIONSHIPS 


From  <A.1> 


<J>  (t,v)  =  SIQ  ®  1 

T,V 

=  f  f  e27TJV'(T~T  )  4> ^ (x-x  1  ,v~v 1 )  Sjq(t',v')  dv*  dx' 


therefore,  if  Eq.  1  of  the  main  text  holds  for  SIQ 

Eri4>0(T,v)|2l  =  f  f  I  /e^^tv,(T'T,)'v,'(T'T“)^E[^(T-T^\)VHI(T*Tl,v■v,)]  | 

L  ^  J  —00  -00  -CO  -OO 


•  E  [ S*0(t" ,v")  Sio(t\v')]  dv"  dx"  dv'  dx’ 
/  /  E  [|4>i(x-t',v-v,)|2]  $io(t',v')  dv'  dx' 


$10  ®  E [1 4> r  (T  >v)  I] 

x,v 


Similarly  from  <A.1> 


4>^{x,a)  =  S°  ®  '  <t>°(x,a) 

°  x,a 


=  /  /  <j>j [a1  (x-x 1 ) ,  £r]  SD  (x  '  ,a' )  da*  dx 
-  -  1  u  10 


-00  0 
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therefore  if  Eq.  1  of  the  main  text  holds  for  $\ 


[|4>0(t ,v)  | 2]  = 


-1/1 1  §*] ♦?[«•(«•)•  b\) 

•  E  |V*(x",a")  Sj0(t  1  ,a‘ )  J  da"  dx"  da1  dx1 
=  /  /  E  ||<|>5  [a'Cx-x1),  g-r]  | 2}  $D(x',a')  da'  dx' 


f  f  f 


-oo  o  -oo 


=  $jq  ®  1  E  |t)>j(x,a)| 

x,a 


Thus  the  fundamental  relationships  of  Fig.  1  of  the  main  text  are  valid. 
Also  from  <A.1> 

Sftr('r>v)  =  Spr  ®  1  $*p(x,v) 


•  f  f  e^v  ^T_T  ^  SAR(x-x',v~v')  SRr(x' ,v‘ )  dv'  dx'  ; 


therefore,  if  Eq.  1  of  the  main  text  holds  for  S.R  and  S, 


e  fs*c<T,v)  s^fT-.v-)]  -777  7 

L  J  -C*3  -00  -00  -oo 

•  e[sab<t’-t">'-v"'>Wt't”’','v">]  e[sbc(t'">v  '")sbc<t”-v")] 


dv1 "dx 1 "dv"dx" 


=  /  /  $flR(x-x",v-v")  $Rr(x",v")  dv"  dx"  <5(x-x')6(v-v') 


=  $BC  $ab(t,v)<5(x-x')6(v-v') 
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and  Eq.  1  of  the  main  text  holds  for  Similarly  from  <A.1> 

sSc(T,0)  =  s“c  ®  '  S^sfT.a) 
t  ,a 

=i  V sab[“'<t-t,)>  I"]sbc<t'-“')  d“'  dT' 

therefore  if  Eq.  1  of  the  main  text  holds  for  S^B  and  sj^ 


00  00  oo  oo 

f  f  f  f  E 

-«  0  *°  0 


*  E[SBC^,"'a",)SBC(T"’v")]  dot'M  dT'"  da"  dT" 

oo  00  _  _ 


®  ' 

T,a 


$AB(t,cx)  5(T-T,)6(a-o') 


and  Eq.  1  of  the  main  text  holds  for  yj^  Thus  the  fundamental  relation¬ 
ships  of  Fig.  2  of  the  main  text  are  valid. 
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APPENDIX  B 

DERIVATION  OF  THE  DOPPLER  AMBIGUITY  FUNCTION  OF  A  UNIFORM  TRAIN  OF 
IDENTICAL  PULSES  AND  SOME  OF  ITS  PROPERTIES 


By  definition  Eq.  lb): 

oo 

YpT(T,a)  =  /  XPT[a(t-x)]x(t)  dt 


Substitution  of  Eq.  15  of  the  main  text  into  this  result  yields  Eq .  14  of 
the  main  text. 


PRECEDING  PAGE  BLANKING!  LWED 

-  •  _  M  -  - - - —  ‘ 
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The  Fourier  transform  of  Eq.  15  of  the  main  text  is  by  definition 


r?(f»a)  = 


=  1  "f  1  j  e'2,,J  -  k  +  1 ^A)  4>(f,a) 

0  k  =  l  I  P 


-Iirj  ~fAT(  ,  n- 1 i  |  -iirj  ^-f(n-|i  |+l-2k)AT)  n 
=  e  a  I  e  e  a  r°(f,a) 

(  n  k=1  >  D 


sinfV(n- 1  i  I  )a]  .  I  >a  - 

sin  Oa]  e^3  -  e-7rja 

bn_ ML  b-n+l 1 1 

b  -  b'*1 

=  bn"  1 1 1 _1  +  bn"ii|_3  +  ...  +  b^M*3  +  b‘n+|il+1 

-  eTt.j(n-|i  | -1  )a  +  eTrj(n-|1  |-3)a  +  t  +  eirj(-n+|1  |+3)a 

+  e*rrj(-n+|i  |+1  )a 


Note  that  this  is  obviously  a  sum  of  cosines;  furthermore,  comparison  of 

this  identity  with  the  above  equation  for  r!?(f,a)  yields  Eq.  17  of  the 
n.  1  1 
main  text  if  a  =  —  fAT  . 

a 

Now  consider  conditions  that  will  ensure  that  the  y?  will  overlap  by  only 
a  negligible  amount  in  the  ambiguity  function.  Consider  first  i  >  0  and 
a  >  1.  When  the  spread  of  Yp  in  t  is  small  compared  with  AT  , 
y^_ .j  [t - ( i - 1 ) AT ,aj  will  not  overlap  y!?(T-iAT,a)  if 
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for  all  k  .  Use  of  the  smallest  k  on  the  left  side  of  this  inequality 
and  the  largest  on  the  right  side  implies  that  the  condition  for  non-overlap 
is 


On  the  other  hand,  for  i  >  0  and  a  <1  ,  the  largest  k  must  be  used  on 
the  left  side  and  the  smallest  on  the  right  side;  hence  the  condition  is 


These  two  conditions  can  be  combined  to  yield  the  condition 
| ot-1  | /a  <  1  / ( n-1 ) .  For  i  <  0  the  same  analysis  leads  to  the  constraint 

| ot-1 1 /a.  <  1/n.  Since  the  latter  constraint  is  more  stringent  than  the 

former,  it  is  the  applicable  constraint. 
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APPENDIX  C 

DERIVATION  OF  THE  APPROXIMATE  MATCHED  FILTER  BANKS 


From  Eq.  9  of  the  main  text  for  xp(t)  =  Xpf^t) 

<t>PT(x,v)  =  —  E  /  e"27TJV(t_T)  X*[t-x-(i-— )AT]  y(t)  dt. 
n  v^T  P  2 

But 

XFM(t)  -n-  (t;T/2)  e+2^fot+*  kt^  , 

implies  Eq.  22  of  the  main  text;  hence 

_i  vf 

1  n  00  -27rj[v(i  £— )  AT  -  — p—  ]  * 

*PT(t.v)  3  —  iJi  fj  XFM[t+v.T.(i.ntL)4T]  y(t)  dt 

Substitution  of  the  definition  Eq.  24  of  the  main  text  of  4>pM(x,o)  into 
this  result  yields  Eq.  23  of  the  main  text. 

Also  from  Eq.  9  of  the  main  text  for  xp(t)  =  Xp^( t ) 

n  co 

^pT(x,a)  =  ~  Z  f  xD  [a(t-T)  -  (1-— )AT]  y(t)  dt  . 


But  the  definition  of  XpM(t) 
hence  1 


n  i  n  «  +2ttj 

<t»pT( t ,a)  =  -  Z  f  e 

/”n  i=l 


given  above  implies  Eq.  25  of  the  main  text; 

(^o_ 

XpM[t+(a-l)f0/k  -x-(i-^)AT/a]y(t)  dt  . 


Substitution  of  the  definition  Eq.  24  of  the  main  text  of  <j>p^(x,o)  into 
this  result  yields  Eq.  26  of  the  main  text. 
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COMMARAIRMED 
C INCHAN 
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